Parametric correlations versus fidelity decay: the symmetry breaking case 
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We provide formulae for fidelity decay and parametric energy correlations for random matrix 
ensembles where time-reversal invariance of the original Hamiltonian is broken by the perturbation. 
Like in the case of a symmetry conserving perturbation a simple relation between both quantities 
exists. Fidelity freeze is observed for systems with even and odd spin. 
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I. INTRODUCTION 

Fidelity presently attracts considerable attention in di- 
verse fields like quantum information, quantum chaotic 
systems and others [l|, [I]- It measures the change of 
quantum dynamics of a state under a modification of the 
Hamiltonian. In quantum information, fidelity measures 
the deviation between a mathematical algorithm and its 
physical implementation. 

Since fidelity requires knowledge of the entire wave 
function for the original and for the modified system 
a measurement of fidelity is a notoriously difficult task. 
However a number of experimental results have been ob- 
tained in microwave billiards, where the perturbation was 
achieved by varying some geometric parameter. There 
are two qualitatively different ways to do this, either by 
a global perturbation, e.g. by moving one wall s'], or 
a local perturbation, e.g. by varying the position of an 
impurity For the first case random matrix theory 
is applicable, and indeed a perfect agreement between 
experiment and theory was found Q. 

On the other hand statistical properties of energy cor- 
relations between spectra of complex quantum systems 
which differ by a parameter-dependent variation have 
been studied experimentally and theoretically Q. This 
quantity can be obtained with great accuracy from scat- 
tering experiments by analyzing the fluctuations of the 
resonances in the scattering cross-section 0| ■ 

From an experimental point of view it is interesting 
to relate fidelity with spectral quantities. This allows an 
indirect measurement of fidelity via an analysis of the 
(parametric) scattering data and the problem of measur- 
ing the entire wave function is circumvented. 



A simple differential relation between fidelity decay 
and parametric energy correlations was established in 
the case that the parameter dependent perturbation falls 
into the same symmetry class as the unperturbed system 

Q . This differential relation was derived earlier in en- 
ergy space by Taniguchi, Simons and coworkers [lol [llj | 
and it was identified with a continuity equation of the 
Calogero-Moser-Sutherland model [12] ■ In Ref. [l^ sim- 
ilar expression were derived for parametric energy cor- 
relations in the case where the perturbation breaks the 
global symmetry of the orgininal unperturbed system. 

Recently billiard experiments could be performed in 
microwave resonators, where time reversal symmetry 
(TRS) was broken by a piece of ferrite 1J| which plays 
the role of the perturbation. From the experimental re- 
sults 5'-matrix elements could be determined and an esti- 
mate of the strength of the TRS-breaking was made. The 
experimental setup seems adequate for a measurement of 
parametric energy correlations and of fidelity decay by a 
TRS breaking perturbation. 

In this paper we therefore analyze the expressions 
found in Refs [l^, [II] for TRS-breaking perturbations 
under the aspect of fidelity and provide formulae for fi- 
delity and parametric form factor as well as differential 
relations between them and discuss their consequences. 



II. DEFINITIONS AND RESULTS 

Fidelity amplitude is defined as a functional of the ini- 
tial wave function. In an ergodic situation it seems rea- 
sonable to replace a specific initial state by a random one. 
In Ref.[lH] the corresponding random matrix model for 
the fidelity amplitude was defined by (?i = 1) 
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The Fourier transform of parametric energy correlations 
is defined by 

X(A[|,Ai,0 = ^(trexp(iii?)trexp(-ztFo)) . (2) 

It was named cross form-factor in Q. The brackets de- 
note an ensemble average. The perturbed Hamiltonian 
H is given as 



H = iJo + A||V|| +AiKL. 



(3) 



Let us first discuss the unperturbed Hamiltonian. We 
assume that for the unperturbed system Hq TRS is con- 
served. The time reversal operator T acts differently 
on systems with integer spin and on systems with half- 
integer spin [l6|. For even spin 7i = C, where C is the 
complex conjugation operator. In this case (called case I 
in the following) Hq is chosen from the ensemble of real 
symmetric matrices, called the Gaussian orthogonal en- 
semble (GOE, 13 — 1). For odd spin systems T2 — CJ, 
where J acts via conjugation with the symplectic met- 
ric. In this case (case II) Hq is chosen from the Gaus- 
sian symplectic ensemble (GSE, /? = 4), consisting of all 
Hermitean matrices which are invariant under T2. The 
ensembles are defined by the averages 



N 

^ {5ii5jk + Sik^ji) GOE 
— {6ii5jk - ^SikSji) GSE . 



.(4) 

For the GSE the matrix entries are quaternions. Zirn- 
bauer and Altland classified random matrix ensembles 
along Cartan's classification of symmetric spaces with 
curvature zero [13, Ell- They called the GOE of Al-type 
and the GSE of All-t^^e. 

In contrast to Ref. [l^ we now assume that the pertur- 
bation contains two parts. One part, named Vji , shares 
the symmetry of Hq and is taken either from a GOE or 
a GSE 



{iv\\hjm)ki) = 



SiiSjk + SikSji GOE 
SiiSjk 



^SikSji GSE . 



(5) 



The second part, V±, breaks the symmetry of Hq and is 
taken for case I from the Gaussian ensemble of matrices 
which change sign under time reversal 7i, being anti- 
symmetric matrices with purely imaginary entries. The 
ensemble is of type B in the classification of For 
case II it is taken from the ensemble of matrices which 
change sign under T2. They are of the block form 



' A 


B ' 




-A* 



a^aK b = b'^ 



(6) 



where A and B are N/2 x N/2 matrices {N even). The 
corresponding ensemble is termed C-type in [isj . The B- 
type and C-type ensembles are defined by the averages 



{{Vi_U{v^)ki) 



SuSjk - kkSji B-type 
5ii6jk + \6ik5ji C-type . 



(7) 



The variance of the matrix elements has been chosen to 
have a mean level spacing one, Z? = 1, for Hq and to be 
of order 1///V for V\\ and V^. 

One parameter A|| characterizes the strength of the per- 
turbation, which conserves TRS. The second one, A_l, is 
the strength of the TRS breaking perturbation. Thereby 
we consider a much wider class of TRS breaking Hamil- 
tonians as before. Observe that for A|| = Aj_ this corre- 
sponds to a perturbation by a Hermitean matrix, i. e. to 
a perturbation which is taken from the Gaussian uni- 
tary ensemble (CUE). This ensemble is called type A 
in [I8]. Thus time reversal symmetry breaking can oc- 
cur in different ways. Symbolically we may write the 
1. h. s. of equation (jS]) as AI -I- A|| AI + A_lB (case I) or as 
AH -I- A|| AI -I- AxC (case II). Usually only the transition 
AI -1- AA is considered, when time-reversal invariance is 
discussed [20| . 

Analysing equation (4) and the following ones of 
Ref. UJ, we find expressions for fidelity amplitude and 
for cross-form factor. To present them concisely we de- 
fine for case I the function 

r u 

zW(A||,Ax,r) j dujvdv 

Max(0,T-l) 

l + Air^Xlir'^ -v^) {t-u){1-t + u) 

2\2 



, (8) 



and for the case II the function 
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where t is time measured in units of Heisenberg time 
tn ~ 2-k/D. Then in the large TV-limit the fidelity as 
defined in Eq. ([!]) is given in both cases by 
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The cross form-factor is given by 



-Z(A,i,A^,r) . (10) 



if(A||,Ai,r) = -r2z(A||,Ai,r) . (II) 

From this follows the relation between fidelity and cross 
form-factor lOl 



/(A||, A^,t) 
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K{\,\ut) . (12) 



This relation can be derived through an universality 
argument without going through a lengthy supersym- 
metric calculation and comparing results. In appendix 
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FIG. 1: Ensemble average of the fidelity amplitude /(A, A,r) 
(solid line) with Ho taken from the GOE (case I) and the 
perturbation taken from the GUE for different perturbation 
strengths A. The results /(A, 0, r) (a pure GOE perturbation, 
dashed lines) and /(O, A, r) (a purely perpendicular pertur- 
bation, dashed-dotted lines) are shown as well for the same 
parameter. 



El we present this derivation extending the method of 
Refs. [1, 1^ to the case of TRS breaking. 

Some details on the the derivation of eqiiations ([5]) to 
(IT2]) from the pertinent formulae of Ref. are given in 
appendix |B] 



III. DISCUSSION 

The double integrals © and © can be evaluated nu- 
merically (see appendix B of Ref. [IJ for a convenient 
parametrization) . Figure [1] shows the fidelity decay in 
case I for different perturbation strengths and for a per- 
turbation taken from a GUE (AI + AA), from a GOE 
(AI-|-A|| AI) and from the B-type ensemble of purely anti- 
symmetric matrices (AI -I- Aj^B). For weak perturbations 
there is nearly no difference between the fidelity decay 
with a GOE and a GUE perturbation, respectively. This 
feature is closely related with the extremely weak fidelity 
decay for an imaginary antisymmetric perturbation. The 
latter was called fidelity freeze and was discussed in 
the context with random matrix theory in Ref. [l^ . The 
diagonal elements of the perturbation in the eigenbasis of 
the original Hamiltonian cause a Gaussian decay, which 
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FIG. 2: Fidelity amplitude /(Ay , Ax, t) for case I (upper 
picture) and for case II (bottom picture) for different values 

of A|| and \i_ and for fixed overall perturbation A = ^/-^y + ^\ 
= 0.1. The values of Ay and Ax are given by the five ratios 
Aj|/Ai = oo (full line), 2 (dotted line), 1 (dashed line), 1/2 
(dashed-dotted line), (dashed-dotted-dotted-dotted line). 



dominates for times larger than Heisenberg time. It was 
therefore predicted [l^ that fidelity decay is much slower 
for perturbations which are purely off-diagonal in the 
eigenbasis of the original Hamiltonian. 

With increasing perturbation the decays for the GOE 
and the GUE perturbation separate, and the freeze be- 
havior get lost. For strong perturbations a recovery of fi- 
delity at Heisenberg time is seen. This is already known 
from [Hi where the cases A + AA and AI -I- AAI were 
discussed. 

For small perturbations and for times much smaller 
than Heisenberg time fidelity decay is governed by 
Fermi's golden rule. In this regime the crucial param- 
eter is A^ = A| + A'^ which is related to the spreading 

width r = 2-k\^D of an unperturbed state. This result 
holds independently of the universality class of the back- 
ground. It is therefore interesting to look on the fidelity 
amplitude for fixed A but different ratios between orthog- 
onal and parallel perturbation. 

In figure [2] fidelity amplitude is plotted for small per- 
turbation strength A = 0.1 and for different ratios be- 
tween A II and Aj^ for case I and case II. 

In case I we see that fidelity amplitude is a monotonous 
function of this ratio for all times. The slowest decay 
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happens for A|| = 0, i. e. when the perturbation in the 
direction of Hq is zero (freeze). In case II fidehty shows 
quahtatively the same behavior, i. e. a slower decay for 
perpendicular perturbations for times beyond Heisenberg 
time. This suggest to define for a general perturbation 
V 



tiHoV 







(13) 



as condition for a fidelity freeze, which is slightly more 
general than the one proposed in However in 

case II the freeze is much less pronounced than in case 
I, indicating that the diagonal elements of V±, albeit 
trVj^ = tiV^Ho — 0, have some impact on the decay. 

A careful look reveals that for times beyond Fermi's 
golden rule but smaller than Heisenberg time in case II 
fidelity decay is slower for a parallel perturbations than 
for a perpendicular perturbation. 

This becomes evident for strong perturbations. In fig- 
ure [3] fidelity amplitude is plotted for the same ratios of 
A|| and Aj_ as before but for strong overall perturbation 
A = 1.5. Case I fidelity decay shows monotonous behav- 
ior as a function A||/Ax and fidelity decay is for all times 
slowest for a perpendicular perturbation. However case 
II fidelity decay is more complicated. For times smaller 
than Heisenberg time decay is slowest for a purely par- 
allel perturbation and fastest for a purely perpendicular 
one. At Heisenberg time a pronounced revival is seen 
for a purely parallel perturbation. The peak decreases 
as the share of the perpendicular perturbation increases. 
Finally for a purely perpendicular perturbation there is 
a minimum at Heisenberg time and no revival at all. 

After Heisenberg time things change. Now decay be- 
comes fastest for a purely parallel perturbation with only 
a tiny second revival at twice the Heisenberg time. For 
a purely perpendicular perturbation the freeze behavior 
comes in and at twice the Heisenberg time a sizable re- 
vival occurs, such that just as in case I for long times 
decay is slowest for a purely perpendicular perturba- 
tion. Somewhere between Heisenberg time and twice the 
Heisenberg time the two curves cross. 

To understand this behavior qualitatively, we recall 
two peculiarities of the GSE: first the spectral rigidity 
is much higher than for the GUE or the GOE. It has 
been argued [l^l that the revival at Heisenberg time is 
a signature of the high spectral rigidity. More generally 
high spectral rigidity favors a slow decay. Second the 
eigenvalues of the GSE are two-fold degenerate (Kramers 
degeneracy) . 

Thus a perpendicular perturbation has two effects: 
first it breaks time reversal invariance and drives the GSE 
into a GUE. Since the latter has lower spectral rigidity, 
this has as consequence that the peak at Heisenberg time 
becomes less and lesser pronounced and for times smaller 
than Heisenberg time decay is enhanced by the perpen- 
dicular perturbation. Second it breaks Kramers degener- 
acy, thus the number of independent levels and therefore 
level density and Heisenberg time double. This leads to 
the pronounced peak at twice the (original) Heisenberg 
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FIG. 3: Same as figure [2] but for a strong overall pertur- 
bation A = 1.5. Fidelity amplitude /(Ay , Ax, t) is shown for 
case I (upper panel) and case II (lower panel) and for the ra- 
tios Ay /Ax = oo (full line), 2 (dotted line), 1 (dashed line), 
1/2 (dashed-dotted line), (dashed-dotted-dotted-dotted 
line). In the lower panel for comparison fidelity amplitude 
/(\/2A, r/2) for a GUE with a GUE perturbation is shown as 
well (thinner full line). 



time. A comparison with the plot of fidelity amplitude 
/(V2A, t/2) of a GUE, with a GUE perturbation(A-FAA) 
shows indeed good agreement. 



In figurelHthe cross form factor is plotted in both cases 
for the same five ratios between A|| and Ax as before. 
Qualitatively the behavior is similar to fidelity amplitude. 
In case I the form factor is smallest for a purely parallel 
perturbation for all times. In case II before Heisenberg 
time the form factor is smallest for a purely perpendicular 
perturbation and largest for a purely parallel one. After 
Heisenberg time the order is inverted. At Heisenberg 
time a logarithmic singularity occurs, which is typical 
for the GSE. For strong perturbations the cross form- 
factor develops peaks at Heisenberg time and for case II 
at twice the Heisenberg time (not shown here). It has 
its cause in the algebraic decay of the cross form-factor 
at these specific times Q. At all other times it decays 
exponentially. 
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FIG. 4: Cross form-factor K{X^^, \±,t) for case I (upper 
panel) and case II (lower panel) for different values of Ay and 
X± and for small overall perturbation A = 0.1. The values of 
A|| and Ax are given by the the ratios Ay/A^ = oo (full line), 
2 (dotted line), 1 (dashed line), 1/2 (dashed-dotted line), 
(dashed-dotted-dotted-dotted line) . 



IV. CONCLUSION 

In conclusion we presented the analytic formulae for 
fidelity amplitude and cross-form factor for parametric 
RMT ensembles, where the time reversal invariance of the 
unperturbed system is broken by the perturbation. The 
general perturbation is split into a parallel component, 
sharing the symmetries of the original Hamiltonian and 
a perpendicular component which maximally breaks this 
symmetry. 

Both possibilities of TRS breaking, even spin GOE— ^ 
GUE and odd spin GSE— > GUE, were discussed on equal 
footing. In the first case a strong freeze effect occurs for 
a purely perpendicular perturbation. It can be explained 
by the absence of diagonal elements of the perturbation 
in the eigenbasis of the unperturbed Hamiltonian. In case 
II long time decay is slowest for a purely perpendicular 
perturbation as well. This leads us to propose trHoV = 
as a more general condition for a reduced fidelity decay. 
However in case II the perturbation has diagonal entries 
in the eigenbasis of Hq and the attenuation of decay is 
much less pronounced than in case I. We are reluctant to 
call this behavior "freeze". We propose to call it "weak 
fidelity freeze" . 

The full Hilbert space is involved in the condition 



trHoV = 0. Therefore it is only applicable to fidelity de- 
cay with respect to a random initial state as considered 
here, to which all states of the Hilbert space contribute. 
For an arbitrary initial state this condition will in general 
not suffice to attenuate fidelity decay. 

In the differential relation between fidelity and cross- 
form factor only the parallel perturbation strength en- 
ters. The relation might be verified experimentally for 
instance in a billiard experiment as described in [ij]. It 
might be used to measure fidelity indirectly via spectral 
correlations. 
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Appendix A: Derivation of Eq. (|12p based on a 
universality argument 

In this section we demonstrate on the example of the 
transition GOE — > GUE (case I) how the method of 
Refs. [1, Q can be extended to the case of symmetry 
breaking. We introduce new variables 



A|| 

^11 - y 



All 



(Al) 



For Hq we allow for a general probability measure in the 
GOE universality class and denote it by di'{Ho), while 
those of V|| and V± are Gaussian measures as before (dV|| 
and dV± include the normalization constants). Since the 
probability measure of Hq is assumed to be general, it 
should be typical and free from any special constraint 
besides the matrix symmetry. 
Now we define 



^ a/3, 7(5 



1 



1 



1 \ 



^zi-Ho + XiiV\iJ 



1 



Z2 - Hq - All Vii - iXi_Vi_ 



(A2) 



Introducing delta-distributions of matrix arguments we 
can express Tap.-yS as 



dHidH2 SiHi -Ho + X\iV\i)T^^^^si^3) 



X SiH^"'^ ~Ho- A|| Vj|) 5{H^^'> - A_lKl) 



(I) 
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where iJi is an x real symmetric matrix, H2 is an 
N X N hermitiai 
Imi?2- Moreover 



N X N hermitian matrix and h!^^^ = Rei/2, -^2^'' 



*^ a/3, 7(5 



1 



1 



(A4) 



All three delta-distributions can be Fourier transformed. 
We find 

To^PM = J dAidA2dA3dHidH2 e2-'fAi(ffi-^o+A||V^|) 

xe^^''''^'^"^"'^^''^^ J'^p^^s ■ (A5) 

Here Ai.2,3 are matrices which have the same symmetry 
as their real space counterparts, namely Hi, i?2^'' and 

H2^\ This means Ai and A2 are N x N real symmetric 
matrices and A3 is an A'^ x iV real antisymmetric matrix. 
The integration domain is the real axis for all indepen- 
dent entries of A„, n — 1,2, 3. The expectation value of 
J'ap.-yS can be written as 

{J^ai3,-i5) = J dv{Ho)dV\\dV±:Fai3,iS 

^g-(l/4)trV|f' + (l/4)tryi^ 

= j dv{Hi3)dAidA2dA:idHidH2Ta[j,^s 

^^2irrtr{Ai(//i-_ffo)+A2(-H"^"'--H"o)+A3-H"^^'} 

Here the brackets (. . .) do not simply mean the expecta- 
tion value. Rather (Tap^-ys) is defined to be the expecta- 
tion value of Tai3,-fS- 

Now we introduce the notation 



tr 



92 



N 

E 



92 



92 



(A7) 



Then it follows from partial integrations that 

/ Q2 \ 



tr 



^aa,l3p ) = -(27r) (tr(AiA2)J"Qa,,3/3) 



and 



(A8) 



{J'cc.,fsp) = -(2^)' (tr(Ai - A2fJ'c.a,pf!) . (A9) 



Here repeated indices stand for summations from 1 to N. 

Let us note that a simultaneous shift of Hi and h!^^^ 
in (jA6p induces a shift of Hq. Although such a shift mod- 
ifies the measure di^lHo), the universality of the spectral 



correlation function implies that {J^aa,pp) is asymptoti- 
cally invariant in the limit N ^ 00. Therefore we obtain 
the following estimate 

/ f d , d \ \_ 



tr 



dHi OH 



-(2^)2(tr(Ai+A2)2j-„„,^^) 
From this it follows that 



O.(AIO) 



d 



9(A) 



tr 



52 



OHidH^^'' 



J^aa,Pt3 



= -7r2(tr(Ai -f A2)2j-„„,^^) 

« 0. (All) 

In order to show that the estimate (jAlOp is indeed cor- 
rect, let us pay attention to Eq. (jA6p . Proper unfolding 
of the energy level correlations requires an 0(1) scaling of 
the eigenvalue density of Hq. Each element of the pertur- 
bation is set to be 0(1), because it should equally scale 
as the mean level spacing. When the eigenvalue density 
is scaled as 0(1), since there are N eigenvalues, each 
eigenvalue Eoj of Ho should typically be 0{N). Then 
the RHS of the identity 



tiiHo) 



N 

E 



(A12) 



becomes 0{N^). In the LHS, on the other hand, we have 
0{N'^) terms, each of which is the square of an element 
of Hq. Therefore each element of Hq is estimated as 
0{N^^^). Then the main contribution to the integral 
over the matrix Hq with respect to the measure di/{Ho) 
in equation (|A6[) comes from a region where the elements 
of Ai + A2 are of order 0{N~^/^). Only in that region a 
rapid oscillation of the exponential factor is avoided. 

It can be seen from the Gaussian factor in Eq. (|A6[) 
that the elements of Ai — A2 are scaled as 0(1). Because 
of the identity 

(Ai - A2)2 = -2(AiA2 + A2A1) + (Ai + A2)2, (A13) 

the elements of (Ai — A2)2 are approximated by the ele- 
ments of — 2(AiA2 + A2A1). Hence we find an estimate 

tr(Ai - A2)2 « -4tr(AiA2), (A14) 

which implies Eg. dAlip . We notice that this estimate 
can only be fulfilled when tr(AiA2) is negative. 
On the other hand, we can readily find 



tr- 



92 



dHidH 



tr 



zi-Hi 



Z2 - H2 



and 
92 



(A15) 



dzidz; 



52 



1 



dzidz2 \zi — H 
1 



1 



l/a/3 V^2 - H2/ 
2 / 1 \ 2 

, (A16) 



tr 



zi-Hi 



1 



Z2 - H2 
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so that 



'92 \ 52 

tr 1^J^aa,l3fj ) = -5 ^ {^afi,l3a) ■ (Al7) 



Comparing (jAllI) and (|A17p . we arrive at 



d 



5(A) 



dzidz2 



(^a/3,^a) , (A18) 



which gives the required relation (I12p between the fi- 
dehty and parametric spectral correlation, respectively 
cross form-factor. 



W = 



(A? + A2 + A2 - 2A1A2 - 1)' 



1 + (a2 + a2 + a2 - 2AA1A2 - 1) (B3) 



Here the plus sign applies to case I and the minus sign to 
case II. The parameter n has the value n = 1 (case I) and 
K — 2 (case II). This factor does not appear in THSA, 
however it does appear in Ref. [l^l- We introduced it, 
such that K{t) is related to K{uj) in both cases via 



K{t) = / duje 



(B4) 



Appendix B: Derivation of equations ([8|l and (|12p 
from Ref. ^ 

In Ref. [11], called THSA in the following, the Fourier- 
transform of the cross form-factor was derived as a three- 
fold integral 

K{x,Xo,Xu,uj) = J d\d\id\2We''^ , (Bl) 

where the integration domains are in case I defined by 
A S [—1, 1], Ai G [1,00], A2 e [1,00] and in case II by A S 
[1,00], Ai e [—1,1], A2 e [0,1]. Setting the parameter 
X = Xu/'i, the expressions for F and W (equations (5) 
and (6) of THSA) are given by 



F± = ±Ki7raj(AiA2 — A) ± 



(A? + A2 - A2 - 1) 



±^(2A2A2-A2-A2-A2 + 1) 



(B2) 



In THSA the function W differs in case I and case II by 
a relative minus sign between two summands in the last 
line of equation (|B3p . This seems to be wrong. Moreover 
in the same line the factor 1/k in the second summand 
is missing in THSA. 

Fourier transformation yields (S(r — A1A2/2 + A) in case 
I and (5(r — A + A1A2) in case II, which allows to integrate 
over A. Equations dS]) to are obtained through the 
transformations 



5 (A1A2 



1) 



V = iVA?A2-Ai-A2 + l 

U = A1A2 



V = VAf Ai - Ai - A2 




The parameters are identified as A|| ~ Xo/2 and Aj^ = 
x„/V2. 
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